rcttinn 


AD-A239  248 


U.  Klf 


24.  SECURITY  ClAS 


2b.  01  CLAS  SIP  CATION/ 


IT  DOCUMENTATION  PAGE 


IVI  MARKINGS 


J  DISTRIBUTION  /  AVAILABILITY  OP  REPORT 
APPROVED  FOR  PUBLIC  RELEASE: 
DISTRIBUTION  UNLIMITED 


4.  PERFORMING  ORGANIZATION  R I  PORT  NUMBCR(S) 

ONR  TECHNICAL  REPORT  // 


).  MONITORING  ORGANIZATION  RIPORT  NUMBIR(S) 


U.  NAME  OP  PERFORMING  ORGANIZATION 

UNIVERSITY  OF  UTAH 


6b.  OPPCE  SYMiOL  I  NAM*  OP  MONITORING- ORGANIZATION 


(If  400JN4M4) 


6c  AOORISS  (Off.  Sum.  end  ZVCode) 

DEPARTMENT  OF  CHEMISTRY 
UNIVERSITY  OF  UTAH 


OFFICE  OF  NAVAL  RESEARCH 
CHEMISTRY  PROG 


7b.  AOORf SS  (Oft  SUm,  4 Wtd  TlPCdbl) 


lb  E;i ;  1  wA 


u.  INU6MI  OP  FUNDING  /  SPONSORING 

lb  OMKI  SYMtOL 

ORGANIZATION 

(It  epodceUe) 

OFFICE  OF  NAVAL  RESEARCH 

ONR 

\ 

Be  AOORESS  (Oty,  Sum.  end  ilk  Code) 

800  NORTH  QUINCY  STREET 
ARLINGTON,  VA  22217-5000 

N00014-89-J-1497 


10.  SOU  PCI  OP  PUNOING  NUMMRS 


MOCftAM 

PHOJfCT 

TASK 

ELI  ME  NT  NO. 

NO. 

NO. 

li.  TITLE  (Inchjde  Security  aumtke ootd 

An  Initio  Studies  of  Ground  and  Excited  Electronic  States  of  MgAr,  CdAr,  and  F.eAr 


TT.  1 


12.  PERSONAL  AUTHOR(S) 

Tprrv  A.  Boatz,  Kold  Lars  Bak,  and  Jack  Simons _ 


134.  TYPE  OP  REPORT  1 3b.  TIM*  COVERED  114.  OATE  OP  RIPORT  (Yeer.  Month.  Oey)  |lS.  PAGECOUNT 

TECHNICAL  imu  rn  Julv  26.  1991  I  34 


16.  SUPPLEMENTARY  NOTATION 


Julv  26,  1991 


OSATI  COOES _ 


CROUP  SUB-GROUP 


18.  SUBJECT  terms  (Cononu*  on  worm  it  netetury  4 nd  identity  oy  0 lot*  number) 


19.  ABSTRACT  (Continue  on  reverie  it  neceaery  4 nd  identity  by  block  number) 

(See  reverse.) 


20.  DISTRIBUTION /AVAILABIUTY  OP  ABSTRACT 
□  UNClASSIPHOAINUMtnO  □  SAMI  AS  RPT 


B  NAM* 

JACK  SIMONS 


one  USIRS 


21.  ABSTRACT  SECURITY  CLASUPICAnON 
SIFTED 


(801)  581-8023 


aiTj-iitiwa ,  jfci^iin  ~y  t>"  ■  vil  r. 


ABSTRACT 

^  The  ground  state  (X  *£”)  and  several  excited  state,  (A  3II,  c  C  'll,  D  'S  , 

-  '**-*-..  -  — - -- 

and  E  JT’’)  ^potential  energy  surfaces  for  the  diatomic  molecules  MgAr,  CdAr,  and 
EeAr  have  been  computed  using  complete  active  space  self -consistent  field 
.CASSCF)  wavefunctions  and  valence  double-  and  triple-zeta  quality  basis  sets 
augmented  with  polarisation  and  diffuse  functions.  Pump -and- probe  laser 
experiments  have  examined  the  quenching  of  excited  singlet  states  of  metal-rare 
gas  complexes  such  as  CdXe  to  produce  triplets  that  dissociate  to  ->  metal  ; 
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-j  vf  ntci 


+  ) 


attractive  C  and  repulsive  c /'rT'curves  of  MgAr  and  CdAr  do  not  intersect  in 
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the  energetically  accessible  region  of  the  C\lI$  surface,  unlike  the 
corresponding  curves  for  the  CdXe  diatom.  These  data  are  consistent  with  the 
absence  ofJPj.Mg  and  ^Pj,  Cd  atoms  in  the  MgAr  and  CdAr  experiments, 
respectively.  However,  an  alternative  quenching  mechanism  involving  "vibronic" 

coupling  between  the  C  'n  vibrational  eigenstates  and  the  continuum  eigenstates 
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of  the  underlying  repulsive  3S*  surface  may  be  operative^  J^his  possibility  is 
examined  and  predicted  to  be  unlikely  for  MgAr  (due  to  small  spin-orbit 
coupling)  and  CdAr  (due  to  unfavorable  vibronic  factors).  BeAr,  which  has  yet 
to  be  probed  experimentally,  is  predicted  to  be  bound  by  770  and  900  cm  x  in  the 
D  'S*  and  E  *2*  states,  respectively,  and  to  display  interesting  potential  curve 
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ABSTRACT 


The  ground  state  (X  £*'•  and  several  excited  state  (A  \T,  c  3Sr,  C  il,  D  XE+ , 
and  E  JX+)  potential  energy  surfaces  for  the  diatomic  molecules  MgAr,  CdAr,  and 
BeAr  have  been  cji.  puted  using  complete  active  space  self-consistent  field 
(CASSCF)  wa/ef unctions  and  valence  double-  and  triple-zeta  quality  basis  sets 
augmented  with  polarization  and  diffuse  functions.  Pump-and-probe  laser 
experiments  have  examined  the  quenching  of  excited  singlet  states  of  metal-rare 
gas  complexes  such  as  CdXe  to  produce  triplets  that  dissociate  to  3Pj  metal 
atoms.  This  quench:  ng,  which  is  ietect.  ed  for  CdXe  buc  not  for  CdAr  or  MgAr,  is 
thought  to  occur  via  a  crossing  or  strong  coupling  of  a  repulsive  triplet  curve 
with  an  attractive  singlet  curve.  The  present  work  indicates  that  the 
attractive  C  1I1  and  repulsive  c  JX+  curves  of  MgAr  and  CdAr  do  not  intersect  in 
the  energetically  accessible  region  of  the  C  1II  surface,  unlike  the 
corresponding  curves  for  the  CdXe  diatom.  These  data  are  consistent  with  the 
absence  of  3Pj  Mg  and  3Pj  Cd  atoms  in  the  MgAr  and  CdAr  experiments, 
respectively.  However,  an  alternative  quenching  mechanism  involving  "vibronic" 
coupling  between  the  C  vibrational  eigenstates  and  the  continuum  eigenstates 
of  the  underlying  repulsive  3 T*  surface  may  be  operative;  this  possibility  is 
examined  and  predicted  to  be  unlikely  for  MgAr  (due  to  small  spin-orbit 
coupling)  and  CdAr  (due  to  unfavorable  vibronic  factors).  BeAr,  which  has  yet 
to  be  probed  experimentally,  is  predicted  to  be  bound  by  770  and  900  cm"1  in  the 
D  1S+  and  E  3S+  states,  respectively,  and  to  display  interesting  potential  curve 
intersections . 
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I.  Introduction 


The  alkaline  earth-rare  gas  and  other  similar  van  der  Waals  diatomic 
systems  have  been  the  focus  of  several  recent  experimental1'7  and  theoretical8"14 
studies,  Of  primary  interest  in  these  works  is  the  characterization  of  both  the 
weakly-bound  van  der  Uaals  ground-state  complex3,6"10, 14  and  low-lying  singlet  and 
triplet  excited  state1"7,10"13  potential  energy  surfaces  of  these  molecules. 

Experimental  studies  of  MgAr3  and  CdAr0  (as  well  as  other  similar 
systems“c,2d)  indicate  that  the  binding  on  the  ground  state  potential  energy 
surface  is  primarily  due  to  a  weak  van  der  Waals  interaction  and  produces  large 
equilibrium  internuclear  distances  (' re  =  4.49  and  4.33  A,  respectively.)  In 
contrast,  the  C  1II  and  A  3n  states  of  MgAr  are  bound  by  over  300  cm"1  and  have 
considerably  smaller  re's  (3.27  and  3.63  A,  respectively.)4  Furthermore,  the 
E  3I+  Rydberg  state  (e.g.,  having  appreciable  (3s)1(4s)1  character  for  MgAr)  is 
bound  by  over  1100  cm"1  (see  Figure  la.)  The  presence  of  minima  on  these 
potential  energy  surfaces  can  be  understood  in  terms  of  an  ion- induced  dipole 
model  in  which  the  rare  gas  atom  is  attracted  to  a  partially  shielded  Mg+  ion 
core.  In  the  II  states,  the  II  orientation  of  the  singly-occupied  p  orbital  on  Mg 
allows  for  a  closer  approach  of  the  rare  gas  atom  than  when  the  p  orbital  is  in 
the  2  orientation. 3  Similarly,  in  the  Rydberg  state,  provided  the  approaching 
rare-gas  atom  penetrates  the  spatially  diffuse  Rydberg  orbital  (which  induces  a 
slight  initial  repulsion  at  long  i niM nuclear  distances),  it  will  be  again  be 
attracted  to  a  partially  shielded  Mg*  ion  core. 

In  the  laser  pump-probe  experiments,2,3,6,7  X  A2+  ground-state  metal-rare  gas 
(i.e.,  M«RG)  molecules  are  excited  to  the  *11  surface  and  subsequent  appearance 
of  3Pj  M  atoms  is  monitored.  Production  of  triplet  M  is  thought  to  occur  via 
crossing  of  the  repulsive  c  32+  state  and  the  initially  populated  C  *n  state. 
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That  is,  the  lower- lying  repulsive  c  32+  state  (which  correlates  at  infinite 
separation  with  3P  V  +  :S  Ar)  intersects  an  energetically  accessible  region  of 
the  bound  C  1n  state.  This  supposition  is  primarily  based  upon  experimental 
studies  of  the  CdXe  diatom,  wiv^re  production  of  the  C  AII  state  leads  to 
efficient  production  of  ^Pj  Cd  atoms  via  a  spin-orbit  coupling  between  the  bound 
C  1II  state  and  the  dissociative  c  3S+  state.7  Indeed,  virtually  no  fluorescence 
from  the  C  1n  state  is  observed  at  all  in  CdXe .  However,  in  analogous 
experiments  on  CdAr,0  just  the  opposite  occurs;  namely,  fluorescence  of  the  C  lII 
state  is  prevalent,  with  no  detectable  formation  of  °Pj  Cd  atoms.  This  is 
rather  puzzling,  since  it  seems  reasonable  to  expect  the  potential  curves  and 
the  strength  of  the  spin-orbit  coupling  of  CdXe  and  CdAr  to  be  quite  similar. 
Likewise,  no  production  of  3Pj  Mg  atoms  is  observed  after  the  formation  of  C  *11 
MgAr , 3  although  this  is  less  surprising  since  the  spin-orbit  coupling  in  Mg  is 
considerably  less  than  that  in  Cd. 

Prior  theoretical  studies  of  MgAr  include  the  all-electron  fourth-order 
Moller- Plesset  calculations9  on  the  ground  state  and  the  perturbative 
calculations^0  of  the  dispersion  energies  of  the  lowest  singlet  and  triplet  2, 

II,  and  A  excited  states  (as  well  as  ths  ground  state.)  Similar  calculations  have 
been  performed  for  CdAr12  and  related  systems.12,13  Apparently,  the  only  prior 
calculation  on  BeAr  has  been  a  rather  approximate  calculation  of  its  ground 
state  potential  energy  surface,  at  the  SCF  level. 1A 

The  purpose  of  the  present  work  is  to  investigate  the  validity  of  the 
physical  pictures  introduced  above  and  to  identify  regions  of  surface 
intersections  and/or  avoided  crossings  for  the  molecules  MgAr,  CdAr,  and  BeAr. 

To  this  end,  qualitatively  correct  potential  energy  surfaces  are  calculated  for 
the  electronic  states  of  interest. 
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II.:  Computational  Methods 


In  order  to  obtain  a  balanced  description  of  all  electronic  states 
considered  here,  complete  active  space  self-consistent  field  (CASSCF) 
wave functions  were  employed.  The  active  space  consisted  of  the  two  3s  electrons 
of  Mg,  distributed  in  all  possible  ways  among  the  3s,  3p,  4s,  and  4p  Mg 
orbitals.  This  space  generates  16,  8,  8,  and  8  configuration  state  functions 
(CSFs)  of  12+,  3S+,  1n,  and  3II  symmetry,  respectively.  It  should  be  noted  that 
not  including  any  of  the  electrons  and  orbitals  on  Ar  in  the  active  space 
excludes  from  the  calculation  those  CSFs  required  to  describe  dispersion 
interactions,  which  are  necessary  for  quantitatively  accurate  calculations  of 
weak  intermolecular  interactions  (such  as  van  der  Waals).  Hence,  these 
calculations  cannot  predict  the  presence  of  shallow  wells  in  the  X  x2+,  C  1n, 

A  3n,  and  E  3S+  states.  However,  the  primary  objective  of  this  study  is  not  a 
quant i ta t ive 1 y  accurate  description  of  well  depths  and  bond  lengths;  these  data 
are  much  more  accurately  known  from  interpretations  of  the  experimental  data. 
Rather,  we  are  seeking  a  qualitatively  correct  picture  of  the  individual 
potential  energy  surfaces  which  will  allow  us  to  identify,  among  other  features, 
probable  regions  of  intersections  and/or  avoided  crossings  between  electronic 
states . 

The  basis  sets  used  in  this  study  are  summarized  in  Table  I.,  For  all  three 
molecules,  the  effective  core  potential  (ECP)  and  the  companion  valence  double- 
zeta  basis  set  of  Stevens,  Basch,  and  Krauss15  were  used  for  Ar.  The  valence 
basis  was  augmented  with  a  single  set  of  diffuse  s  and  p  functions  as  well  as 
two  sets  of  five  d-type  polarization  functions.  (See  Table  I  for  the  values  of 
the  exponents  used.)  In  order  to  test  the  validity  of  using  an  effective  core 
potential  on  Ar,  the  calculations  of  all  the  potential  energy  curves  of  BeAr 
were  repeated  using  an  all-electron  treatment  of  Ar.  Virtually  identical 
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potential  energy  curves  are  obtained  for  all  of  the  electronic  states.  For  the 
all-electron  calculations,  the  Ar  basis  consisted  of  the  McLean- Chandler 
(12s9p/6s5p)  set16  plus  the  same  set  of  diffuse  and  polarization  functions  used 
in  the  ECP  computations. 

For  each  of  the  metal  atoms  considered  in  this  work,  the  basis  set  used  is 
constructed  in  a  similar  manner.  Starting  with  a  basis  set  which  is  at  least  of 
valence  double-zeta  quality,  a  single  set  of  five  d-type  polarization  functions 
are  added.  Next,  in  order  to  have  a  basis  set  which  is  sufficiently  flexible  to 
describe  the  triplet  (ns)([n+l]s)  Rydberg  states  (i.e.,  E  3S+)  ,  two  sets  of 
diffuse  s  and  p  functions  are  added  (see  Table  I  for  the  exponents,)  An 
exception  is  made  for  cadmium,  for  which  only  a  single  set  of  diffuse  s  and  p 
functions  were  added.  For  magnesium,  the  McLean-Chandler  (12s9p/6s5p)  basis 
set16  is  used.  The  6-311G  basis  set  of  Pople  et  al.17  is  used  for  beryHium. 

The  ECP  of  Hay  and  Wadt18  and  the  associated  double-zeta  valence  basis  set  is 
used  for  cadmium.  The  MESSKIT19  suite  of  computer  programs  was  used  in 
performing  all  energy  calculations.  All  calculated  potential  energy  curves  were 
corrected  for  basis  set  superposition  error  using  the  counterpoise  method  of 
Boys  and  Bernard!.20 

III.  Results  and  Discussion 

In  order  to  compare  the  present  calculations  with  experiment,  the  atomic 
energies  of  the  metal  atoms  for  the  states  which  correlate  to  the  X  1Z*f  A  3II, 
c  3E+,  C  D  1S+,  and  E  3£+  states  of  the  M»RG  complexes  were  calculated,  using 

the  same  active  space  and  basis  sets  as  outlined  in  the  previous  section.  The 
results  are  summarized  in  Table  II.  With  the  exception  of  the  Be  !S  ->  !P 
transition  energy  (which  is  about  0.9  eV  too  high),  all  calculated  transition 
energies  agree  to  within  approximately  0.5  eV  of  the  experimental  values. 
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For  ease  of  chemical  interpretation,  Table  III  summarizes  the  relationship 
between  the  spectroscopic  notation  for  each  electronic  state  and  the  dominant 
electron  configuration,  in  terms  of  the  orbitals  on  the  natal  atom. 

The  calculated  potential  energy  curves  for  MgAr,  CdAr.  ;:nd  BeAr  are  shown 
in  Figures  la,  lb,  and  lc,  respectively.  (No  local  minima  are  observed  in  any 
of  the  MgAr  and  CdAr  curves,  in  contrast  to  experimental  findings.  As  mentioned 
earlier,  this  is  to  be  expected  based  upon  the  nature  of  the  present 
calculations  and  is  of  minor  significance  in  regard  to  the  objectives  of  this 
paper. ) 

MgAr:  Figure  la  shows  that  the  X  12+  potential  remains  "flat"  until  the 

internuclear  distance  R  decreases  to  about  4  A,  after  which  it  rises  quite 
sharply.  The  A  3II  and  c  3S+  states,  which  are  degenerate  at  R  «  «,  remain 
nearly  so  until  R  decreases  to  about  5  A.  At  this  point,  the  c  3S+  surface 
begins  to  rise  in  energy  while  the  A  3n  state  remains  level  until  R  shrinks  to 
near  4  A.  Similar  behavior  is  observed  for  the  C  1II  and  D  12+  states  (which  are 
also  degenerate  at  infinite  separation,  in  the  absence  of  spin-orbit  coupling), 
although  the  C  1II  state  is  less  strongly  repulsive  than  its  triplet  counterpart. 
These  observations  are  consistent  with  the  fractional  charge- induced  dipole 
model,3  which  states  that  the  II  orientation  of  the  singly-occupied  p  orbital  on 
Mg  allows  for  a  closer  approach  of  the  argon  atom  than  does  the  2  orientation. 
Finally,  the  E  3S+  Rydberg  state  remains  level  until  R  decreases  to  near  3  A, 
and  thus  behaves  similarly  to  the  n  states. 

Of  particular  interest  in  Figure  la  are  the  relative  locations  of  the  C  1n 
and  c  32+  curves.  Based  on  experimental  studies  of  CdXe,7  it  might  be 
anticipated  that  these  two  curves  intersect  In  an  energetically  accessible 
region  of  the  C  XII  surface;  i.e.,  at  an  energy  below  the  asymptotic  limit  of  XP 
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Mg  +  XS  Ar,  Inspection  of  Figure  la  shows  that  these  curves  intersect  well 
above  this  limit.  In  Figure  2a,  the  C  lII  and  c  3S+  curves  have  been  shifted  so 
that  the  asymptotic  *P  Mg  +  lS  Ar  and  3P  Mg  +  lS  Ar  energies  match  the 
experimental  values.  Even  so,  the  intersection  of  these  two  curves  still  occurs 
far  above  the  asymptotic  XP  Mg  +  lS  Ar  limit.  While  it  is  possible  that  the 
level  of  calculation  here  is  inadequate  to  correctly  describe  the  behavior  of 
the  C  !n  and  c  3S+  states  for  sma1!  internuclear  distances,  it  is  also  possible 
that  a  mechanism  other  than  curve  crossing  exists  whereby  a  spin-orbit  assisted 
transition  between  these  two  states  can  take  place.  Although  it  seems  unlikely 
that  these  curves  actually  intersect  at  an  energy  below  that  of  l?  Mg  +  XS  Ar, 
it  seems  quite  probable  that  these  curves  will  closely  approach  each  other  for 
energies  near  this  energy  limit.  This,  in  turn,  could  allow  for  a  large  overlap 
between  the  vibrational  eigenstates  of  the  bound  C  lII  state  with  the 
iscenergetic  continuum  eigenstates  of  the  repulsive  c  3S+  state.  This 
possibility  will  later  be  explored  in  greater  detail  (vide  infra) . 

CdAr :  The  computed  potential  energy  curves  of  CdAr  are  found  in  Figure  lb  and 
are  quite  similar  to  those  of  MgAr  (see  Figure  la).  As  in  MgAr,  the  ground 
state  potential  starts  to  rise  in  the  vicinity  of  R  -  4  A.  However,  both  the 
c  3E+  and  D  lH+  curves  become  repulsive  at  a  slightly  larger  internuclear 
distance  (near  R  -  5  A)  relative  to  MgAr. 

Inspection  of  Figure  lb  shows  that  the  c  3S+  and  C  XII  curves  intersect  at 
an  energy  well  above  the  *P  Cd  +  lS  Ar  limit.  This  is  still  true  even  after 
these  curves  have  been  shifted  so  that  the  asymptotic  energies  match  those  of 
experiment  (see  Figure  2b.)  Since  the  c  3S+  and  C  1II  states  seem  likely  not  to 
intersect  at  an  energy  below  1 P  Cd  +  1S  Ar,  the  same  vibronic  interaction 
mechanism  outlined  above  may  again  be  the  primary  means  of  a  spin-orbit  assisted 
transition  between  these  two  curves  (vide  infra) . 
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BeAr:  Figure  lc  shows  the  calculated  potential  energy  curves  of  BeAr.  These 


surfaces  exhibit  the  same  general  characteristics  of  those  for  MgAr  and  Cdar; 
namely,  that  the  c  3S+  and  D  XE+  curves  become  repulsive  at  longer  R  than  do 
their  n  counterparts  and  that  an  intersection  of  the  c  3S+  and  C  1n  curves  at  an 
energy  below  1P  Be  +  XS  Ar  seems  improbable.  Unlike  MgAr  and  CdAr,  however,  the 
calculated  potential  energy  curves  for  BeAr  show  the  presence  of  local  minima  in 
the  D  XE+  and  E  3S+  curves,  with  well  depths  of  approximately  770  and  900  cm"1, 
respectively.  The  former  apparently  is  caused  by  an  avoided  crossing  with  a 
higher  state  of  symmetry  (not  shown  in  Figure  lc.)  The  presence  of  a  well  on 
the  E  3Et  Rydberg  surface  is  consistent  with  the  fractional  charge  -  induced 
dipole  model  as  discussed  in  the  Introduction.  The  "shoulder"  on  the  repulsive 
limb  of  the  c  3S+  curve  is  apparently  the  result  of  an  avoided  crossing  with  the 
higher- lying  E  3S+  state.  Note  that  in  the  region  of  this  avoided  crossing,  the 
c  3E+  state  is  diverted  away  from  the  C  1II  curve.  Consequently,  if  the  latter 
two  curves  intersect,  the  point  of  intersection  will  occur  at  a  smaller 
internuclear  distance  (and  hence  at  a  higher  energy)  than  if  the  avoided 
crossing  was  absent. 

As  seen  in  Figure  lc,  the  calculated  c  3S+  and  C  1II  curves  of  BeAr  do  not 
intersect  below  the  XP  Be  +  XS  Ar  energy.  However,  when  these  two  curves  are 
shifted  to  match  the  corresponding  experimental  asymptotic  energies ,  their 
intersection  occurs  at  an  energy  near  the  asymptotic  limit  (see  Figure  2c.)  This 
is  in  contrast  to  the  analogous  curves  of  MgAr  and  CdAr  and  suggests  that,  in 
the  case  of  BeAr,  predissociation  to  the  c  32+  state  from  the  C  1  II  surface  may 
indeed  take  place  via  an  actual  curve  crossing,  rather  than  by  a  vibronic 
mechanism  (vide  infra.) 

In  Figure  2c,  all  of  the  calculated  potential  energy  curves  of  BeAr  have 
been  shifted  so  that  their  asymptotic  energies  match  the  corresponding 
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experimental  atomic  energies  of  the  Be  atom.  We  believe  that  this  gives  a 
reasonably  accurate  picture  of  the  ground  and  low-lying  excited  states  of  BeAr, 
and  so  the  potential  energy  surfaces  in  Figure  2c  will  form  the  basis  of  our 
discussion  on  the  experimental  implications  of  the  wells  on  the  D  1S+  and  E  3Z+ 
surfaces . 

The  presence  of  these  wells  is  intriguing  and  suggests  that  these  excited 
states  of  BeAr  may  be  detectable  in  laser  pump-probe  experiments  similar  to 
those  performed  on  MgAr3,4  and  CdAr.5,6  In  particular,  if  a  metastable  triplet 
BeAr  complex  on  the  (presumably  bound)  A  3II  surface  can  be  prepared  in  a 
supersonic  jet  expansion  (as  is  done  for  MgArA)  ,  then  laser  excitation  to  the 
bound  E  3Z+  Rydber'  state  should  be  possible.  Furthermore,  once  triplet  Rydberg 
BeAr  has  been  prepared,  several  routes  may  then  be  taken  to  other  states.  Since 
the  c  3Z+  and  E  3Z+  states  undergo  an  avoided  crossing  in  the  vicinity  of  the  E 
3S+  minimum,  a  radiationless  transition  back  to  the  c  3Z+  surface  is  possible. 
Note,  however,  that  the  degree  to  which  these  two  states  avoid  each  other  is 
probably  overemphasized  in  Figure  2c.  This  is  because  the  energy  spacing  of  the 
corresponding  unshifted  curves  in  Figure  1c  (from  which  the  curves  in  Figure  2c 
are  obtained)  is  somewhat  smaller  than  indicated  by  the  experimental  Be  atomic 
excitation  energies. 

Because  the  spin-orbit  coupling  is  small  for  beryllium,  a  spin- forbidden 
transition  from  the  E  3Z+  Rydberg  state  of  BeAr  to  the  underlying  D  1Z+  state  is 
unlikely  to  occur  to  a  significant  extent.  However,  such  a  transition  could 
allow  for  the  characterization  of  both  the  short-range  attractive  and  long-range 
repulsive  regions  of  the  D  x2+  surface. 

Vibronic  Coupling:  Although  the  c  3Z+  and  C  xn  curves  of  MgAr  and  CdAr 
apparently  do  not  intersect  at  an  energy  below  the  l?  Mg  +  lS  Ar  and  1 P  Cd  +  XS 
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Ar  asymptotic  limits,  respectively,  there  still  may  exist  other  means  by  which 
these  two  potential  energy  surfaces  can  interact.  Since  these  two  curves  are  in 
close  proximity  at  small  internuclear  distances,  a  possible  mechanism  is  for  a 
"vibronic"  coupling  between  the  bound  vibrational  eigenstates  of  the  C  ln  state 
and  the  continuum  eigenstates  of  the  unbound  c  3Z+  state.  This  kind  of  coupling 
is  illustrated  in  Figure  3,  which  shows  a  vibrational  eigenfunction  of  a  Morse 
potential21  and  the  isoenergetic  continuum  eigenfunction  of  the  underlying 
repulsive  triplet  curve  superimposed  on  their  respective  potential  energy 
surfaces.  The  net  overlap  of  these  two  eigenfunctions  serves  as  an  indication 
of  the  likelihood  of  a  connecting  vibronic  interaction.  Since  the  continuum 
eigenstate  oscillates  much  more  rapidly  than  does  the  Morse  potential 
eigenstate,  the  overlap  of  these  two  functions  is  expected  to  be  virtually  zero, 
with  the  possible  exception  of  the  region  of  the  repulsive  "wall"  where  the  two 
potential  energy  curves  are  in  close  proximity.  In  this  region,  the  rate  of 
oscillation  of  the  continuum  eigenfunction  is  somewhat  reduced,  which  may  allow 
for  a  nonzero  overlap  with  the  innermost  "envelope"  of  the  Morse  potential 
eigenfunction. 

In  order  to  investigate  the  likelihood  for  such  vibronic  couplings  in  MgAr 
and  CdAr,  the  C  1II  states  of  these  diatomics  were  modelled  using  Morse 
potentials,  with  the  Morse  parameters  determined  from  fits  to  experimental 
data.3,6  Our  calculated  triplet  curves  in  Figures  2a  and  2b  were  used  to 
represent  the  underlying  c  3Z+  states.  Then,  for  each  experimentally  observed 
vibrational  level  of  the  bound  C  1n  state,  the  corresponding  Morse  vibrational 
eigenfunction  was  computed.  At  the  corresponding  energy  on  the  underlying  c  3Z+ 
state,  the  continuum  nuclear  wavefunc Lion  was  then  calculated.  (Since  the 
continuum  eigenfunctions  of  MgAr  and  CdAr  do  not  approach  zero  in  the  limit  of 
large  internuclear  distances,  they  were  "box-normalized"  using  an  arbitrary 
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length  L  which  extended  well  beyond  the  outer  envelopes  of  all  of  the  Morse 
potential  eigenfunctions  calculated  here.)  Finally,  the  overlap  was  computed  via 
numerical  integration  of  the  product  of  these  two  nuclear  wavefunctions ♦ 

In  MgAr,  there  are  eight  (v-0-7)  vibrational  levels  observed  in  the  C  1II 
state.3  Figures  4a  and  4b  are  plots  of  the  v-0  and  v-7  vibrational 
eigenfunctions,  respectively,  of  the  C  1II  state  of  MgAr  superimposed  with  the 
corresponding  isoenergetic  continuum  eigenfunctions  of  the  underlying  unbound 
c  3S+  state.  Of  particular  importance  in  these  plots  is  the  relative  position 
of  the  innermost  envelopes  of  the  Morse  eigenfunction  and  the  continuum 
eigenfunction,  since  the  total  overlap  is  determined  primarily  by  the  degree  of 
overlap  of  these  two  envelopes.  For  the  v-0  level  of  C  1II  MgAr  (Figure  4a),  the 
innermost  envelope  of  the  continuum  eigenfunction  occurs  well  before  that  of  the 
Morse  eigenfunction.  Therefore,  a  small  overlap  is  anticipated  for  this 
vibrational  level.  However,  for  v-7  (Figure  4b),  one  observes  a  significant 
non-zero  overlap  of  the  respective  innermost  envelopes  and  therefore  a 
relatively  large  total  overlap  is  expected. 

The  results  of  the  numerical  integrations  of  the  Morse  potential 
eigenfunctions  with  the  corresponding  continuum  eigenfunctions  for  all  eight 
vibrational  levels  of  MgAr  are  summarized  in  Table  IV.  There  are  two  sets  of 
values  listed  in  Table  IV.  One  set  consists  of  the  actual  overlap  magnitudes 
and  the  other  set  is  the  percentages  of  the  maximum  possible  overlaps.  The 
maximum  overlap  is  defined  as  the  overlap  which  occurs  when  the  c  3E+  state 
intersects  the  C  lTl  state  at  the  energy  of  the  particular  C  iII  vibrational 
level .  The  data  in  Table  IV  shows  that  both  the  overlap  magnitudes  and 
percentages  increase  with  increasing  vibrational  quantum  number..  This  implies 
that  the  C  XII  and  c  3S+  curves  are  becoming  closer  as  the  vibrational  energy  of 
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the  C  1II  state  increases  and  that  there  is  a  concomitant  increase  in  the 
likelihood  of  a  vibronic  coupling  of  these  two  states. 

There  are  twelve  observed  vibrational  levels  (v=0-ll)  for  CdAr.6  Figures  5a 
and  5b  are  plots  of  the  v-0  and  v-11  vibrational  eigenfunctions,  respectively, 
of  the  C  1rt  state  superimposed  on  the  corresponding  continuum  eigenfunctions  of 
the  underlying  c  3Z+  state.  Note  that  the  oscillations  of  the  continuum 
eigenfunctions  are  more  rapid  for  CdAr  than  for  MgAr.  This  is  due  to  the  larger 
XP  -  3P  energy  difference  in  Cd  (1.69  eV)  relative  to  Mg  (1.65  eV) .  Similar  to 
MgAr,  the  v-0  plot  (Figure  5a)  shows  virtually  zero  overlap  of  the  innermost 
envelopes  of  the  respective  eigenfunctions.  However,  even  for  v-11,  there  is 
still  very  little  overlap  of  the  inner  envelopes  (see  Figure  5b).  Table  V 
summarizes  the  computed  overlaps  and  the  maximum  overlap  percentages  for  each  of 
the  twelve  observed  vibrational  levels  of  C  XII  CdAr.  As  in  MgAr,  both  the 
overlap  magnitudes  and  percentages  increase  with  increasing  vibrational  quantum 
number.  However,  both  sets  of  data  are  significantly  smaller  than  their  MgAr 
counterparts.  This  indicates  (ignoring  the  difference  in  magnitude  of  Cd  and  Mg 
spin-orbit  couplings)  that  CdAr  is  less  likely  to  undergo  a  vibronic  coupling  of 
the  C  xn  and  c  3S+  states  than  is  MgAr. 

IV.  Conclusions 

This  paper  has  attempted  to  provide  a  detailed  look  at  the  ground  state  and 
several  of  the  excited  state  potential  curves  of  the  diatomic  molecules  MgAr, 
CdAr,  and  BeAr.  In  contrast  to  implications  from  experimental  studies  of  CdXe, 
the  computed  potentials  show  no  intersection  of  the  C  1n  and  c  3S+  states  in 
regions  below  the  energetically  accessible  :P  Mg  (or  1P  Cd)  +  1S  Ar  limit.  A 
mechanism  involving  "vibronic"  coupling  of  the  vibrational  eigenstates  of  the  C 
xn  state  and  the  continuum  eigenstates  of  the  underlying  c  3Z+  state  is  proposed 
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as  an  alternative  means  of  intersystem  crossing  in  MgAr  and  CdAr.  Calculation 
of  the  overlap  of  the  nuclear  wavefunctions  of  the  c  3S+  and  C  1n  states  suggest 
that  this  type  of  vibronic  coupling  is  more  likely  to  take  place  in  MgAr  than 
CdAr;  of  course,  spin-orbit  coupling  is  larger  in  the  CdAr  case  than  for  MgAr. 
The  overall  intersvstem  crossing  rate  is  expected  to  depend  on  the  product  of 
spin-orbit  factors  and  a  vibronic  factor  and  thus  to  be  small  for  both  MgAr  (due 
to  the  small  spin-orbit  factor)  and  CdAr  (due  to  the  small  vibronic  factor)* 

The  D  lE+  and  E  3S+  states  of  BeAr  are  predicted  to  be  bound  by  770  and  900 
cm**1,  respectively,  and  thus  may  constitute  interesting  synthetic  targets  for 
the  experimentalist . 
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Figure  Captions 


la.  Calculated  potential  energy  surfaces  of  MgAr,  with  the  experimentally 
determined  equilibrium  bond  lengths  denoted  "Re"  and  the  ground-state 
v~0  turning  points  shown  by  (), 

lb.  Calculated  potential  energy  surfaces  of  CdAr,  with  the  experimentally 
determined  equilibrium  bond  lengths  denoted  "Re”  and  the  ground-state 
v-0  turning  points  shown  by  (). 

lc.  Calculated  potential  energy  surfaces  of  BeAr. 

2a.  Calculated  C  1II  and  c  3S+  potential  energy  surfaces  of  Mger, 

shifted  to  asymptotically  match  the  experimental  energies  of  :P  and  3P  Mg, 
respectively . 

2b.  Calculated  C  lR  and  c  32+  potential  energy  surfaces  of  CdAr, 

shifted  to  asymptotically  match  the  experimental  energies  of  XP  and  3P  Cd, 
respectively . 

2c.  Calculated  potential  energy  surfaces  of  BeAr,  shifted  to  asymptotically 
match  the  experimental  energies  of  ,  3P,  and  3S  Be. 

3.  "Vibronic”  coupling  between  a  vibrational  eigenstate  of  a  bound 

lII  state  and  the  isoenergetic  continuum  eigenstate  of  an 
underlying  dissociative  3S+  state. 

4a.  Morse  potential  vibrational  eigenfunction  for  the  v«0  vibrational  level 
of  C  1II  MgAr  and  th  isoenergetic  continuum  eigenfunction  of 
the  underlying  c  3I  state, 

4b.  Morse  potential  vibrational  eigenfunction  for  the  v-7  vibrational  level 
of  C  1II  MgAr  and  the  isoenergetic  continuum  eigenfunction  of 
the  underlying  c  3E+  state. 

5a.  Morse  potential  vibrational  eigenfunction  for  the  v~0  vibrational  level 
of  C  1n  CdAr  and  the  isoenergetic  continuum  eigenfunction  of 
the  underlying  c  3S+  state. 

5b.  Morse  potential  vibrational  eigenfunction  for  the  v-11  vibrational  level 
of  C  1II  CdAr  and  the  isoenergetic  continuum  eigenfunction  of 
the  underlying  c  3S+  state. 


Table  I.  Summary  of  Basis  Sets  and  Effective  Core  Potentials. 


Ar 

Ard 

Mg 

Cd 

Be 

ECPa 

SBK 

None 

None 

HW 

None 

Basis'5 

SBK  DZ 

MC 

MC 

HW  DZ 

6  - 311GC 

<*s 

0.027 

0.027 

0.0209, 

0.00825 

0.02176 

0.0285, 

0.01125 

qp 

0.021 

0.021 

0.0152, 

0.0060 

0.0162 

0.0171, 

0.00675 

Qd 

1.70, 

0.425 

1.70, 

0.425 

0.175 

None 

0 . 400 

a  i.SBK.. 

("Htf")  denotes 

the  ECP  from 

reference 

15  (18). 

b  "SBK  DZ"  ("HW  DZ")  denotes  the  double-zeta  valence  basis  defined 

in  reference  15  (18).  "MC"  denotes  the  McLean-Chandler  (12s9p/6s5p)  basis 
from  reference  16. 
c  See  reference  17. 

d  Used  only  in  the  BeAr  all-electron  calculations. 


Table  II.  Calculated  and  .  ..perimental  Atomic  Excitation  Energies.1 


Atom 

lS  ->  3P 

3S  ->  3P 

!S  ->  3S 

Mg 

2.71 

4.78 

4.98 

(2.69) 

(4.34) 

(5.12) 

Cd 

3.38 

5 . 31 

5.85 

(3.73) 

(5.42) 

(6.37) 

Be 

3.08 

6.20 

6.37 

(2.73) 

(5.29) 

(6.46) 

a  In  eV. 

Experimental  values  (taken  from  Moore, 

C.E.  "Atomic  Energy 

Levels 

",  Vols.  1  and  3, 

Nat.  Stand.  Ref.  Data. 

Ser.,  NBS ,  1971)  are 

in  parentheses . 

Table  III.  Dominant  Electron  Configuration  for  Each  Electronic  State. 


Electronic  State  Configuration 


x 

a  3n,  c  :n 

c  3s\  d  1s1 
E  3S+ 


(ns)2 

(ns):(npff): 
(ns)1(np<7)1 
(ns) l( [n+1 Is)1 


Table  IV.  Values  of  Nuclear  Wavefunction  Overlap  Integrals  of  MgAr. 


Vibrational 
Quantum  Number 


| <c  3S+  I  C  MI  >|b  Percent 

of  Maximum  Overlap13 


0 

8 . 3xlO'5 

0.10 

1 

2.6x10"“ 

0.36 

2 

5.8x10"“ 

0.87 

3 

l.lxlO"3 

1.70 

4 

1 . 8xl0‘3 

2.92 

5 

2 . 6xl0'3 

4.52 

6 

3 . 6xl0'3 

6.52 

/ 

4.6xl0"3 

S .  82 

a  All  values  are  unitless. 
b  See  text. 


Table  V.  Values  of  Nuclear  Wavefunction  Overlap  Integrals  of  CdAr.a 


Vibrational 
Quantum  Number 

l<c  3s+  |  c  1n  > | b 

Percent 

of  Maximum  Overlap1 

0 

5 . 9xl0'7 

0.00 

1 

4 . 2xl0"7 

0.00 

2 

1.5xl0‘6 

0.00 

3 

2 . 3xl0"6 

0.00 

4 

3 . 4xl0'6 

0.01 

5 

6 . 7xl0"6 

0.01 

6 

l.lxlO'5 

0.02 

/' 

1 . 8xl0‘5 

0.04 

8 

2 . 6xl0"5 

0.06 

9 

3 . 6xl0"s 

0.08 

10 

1 . 7xlO*5 

0.11 

11 

5 . 9xl0"5 

0.15 

a 


b 


All  values  are  unitless. 
See  text. 
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